Twin Boundaries in <i-wave Superconductors 
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Twin boundaries in orthorhombic d-wave superconductors are investigated numerically using the 
Bogoliubov-deGennes formalism within the context of an extended Hubbard model. The twin 
boundaries are represented by tetragonal regions of variable width, with a reduced chemical po- 
tential. For sufficiently large twin boundary width and change in chemical potential, an induced 
s-wave component may break time-reversal symmetry at a low temperature T* . The temperature 
T* , and the magnitude of the complex component, are found to depend strongly on electron density. 
The results are compared with recent tunneling measurements. 
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In spite of mounting experimental evidence that the 
high-temperature superconductors have an order param- 
eter with d x 2_ y 2 (d-wave) symmetry^ a number of exper- 
iments on both twinned and untwinned YBa2Cu307_5 
(YBCO) suggestpthe presence of an additional s-wave 
order parameter.1313 By symmetry, a small s-wave com- 
ponent always coexists with a predominantly d-wave or- 
der parameter in an orthorhombic superconductor such 
as YBCO, and changes its sign (relative to the d-wave 
component) across a twin boundary!] The experimen- 
tal results can be understood either if the s-wave com- 
ponent of the order parameter breaks time-reversal (T) 
symmetry near the twin-boundary at low temperatures, 
as it can near surfaces,ETc2l or if there is far more of 
one kind of .twin domain (i.e. twin boundaries form in 
groups). EJEj Recent SQUID measurements on vortices 
trapped by twin boundaries in YBCO did not detect the, 
fractional flux that would accompany local T-violation.E3 
In the present work, however, we present strong evidence 
that such a symmetry breaking could indeed occur in the 
vicinity of twin boundaries at low temperatures under 
certain conditions. 

Sigrist et alu have addressed the possibility of T- 
violation near boundaries by considering the Ginzburg- 
Landau (GL) free energy for a homogeneous orthorhom- 
bic d-wave superconductor: 

F s = F n + a d \d\ 2 + a s \s\ 2 + PM A + P 2 \s\ A + f3 3 \s\ 2 \d\ 2 
+ /3 4 (s* V + s 2 d* 2 ) + (3 5 (s*d + sd*), (1) 

where s and d are the s-wave and ci-wave components of 
a superconducting order parameter of the form d + e l6 s, 
with 9 the relative phase between s and d. Only the 
lowest-order orthorhombic term is kept in the GL free 
energy since s -C d for small (3$. Assuming all the coeffi- 
cients are positive (except ad), the /3s term favors 9 = utt 
(n integer), while the /?4 term favors 9 = tt/2. Thus, in 
a tetragonal superconductor there can be a continuous 
transition to a bulk T-violating phase at a tempepature 
T* given by a s {T*) = [2/3 4 (T*) - /3 3 (T*)} d 2 (T*)B The 
coefficients of the GL free energy can be chosen such that 
T* < for a uniform system; the orthorhombic term fur- 
ther discourages a positive T*. Yet the suppression of 



d near an inhomogeneity, or the presence of a tetrago- 
nal region within an orthorhombic system, could induce 
local T-violation for a finite temperature < T < T* . 
Furthermore, the associated complex s-wave component 
could be large relative to \d\, and would vary on a new 
length scale. In contrast, an s-wave component nucle- 
ated solely through spatial /variations of d, such as is 
found in magnetic vorticescjtJ or near impurities ,lj is 
usually small relative to \d\ (unless a s — > 0, which only 
occurs for densities-just above a crossover to bulk s-wave 
superconductivitytZl) and varies on a length scale of the 
d-wave coherence length. 

It is not clear, however, whether GL theory (which 
is strictly valid only near the bulk superconducting 
transition temperature) can reliably describe the low- 
temperature regime associated with T-violation. In the 
present work, twin boundaries in d-wave superconductors 
are investigated numerically using Bogoliubov-deGennes 
theory. We employ an extended Hubbard model which 
gives rise to d-w-aye superconductivity in a restricted pa- 
rameter regime.Ea Despite its simplicity, results obtained 
previously using this modell3'E2l are consistent with those 
obtained within the context, rof a model better repre- 
senting—the high-T c oxides, t3'Eil and with experimental 
results.E3E3 Twin boundaries are modeled as tetragonal 
regions of varying widths and reduced chemical poten- 
tial, in order to approximate the experimental observa- 
tions that twin boundaries are oxygen-deficient (i.e. lo- 
cally antifeppmagnetically insulating) regions, generally 
7-40Awide.Efl 

The Hamiltonian for the extended Hubbard model is: 



(ij)a ia ia 

tt n i<T n ja >, (2) 



{ij) c 



where the sums are over spin and nearest-neighbors on 
the square lattice, is a direction-dependent hopping 
parameter used to model orthorhombicity, fj, is the chem- 
ical potential, /i 7 is a site-dependent impurity potential 
representing the depletion of the carrier density at the 
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twin boundary, and Vq and V\ are on-site and nearest- 
neighbor interactions, respectively (V > denotes at- 
traction). Choosing the unit cell as shown in Fig. [I], we 
can exploit the translational invariance of the Hamilto- 
nian in the (110) direction. With R = x + y parallel and 
r = — x + y perpendicular to the twin direction, we obtain 
the Bogoliubov-deGennes (BdG) equations: 



(i 




( u nt k(r a ) \ 






\v n ,-k(r a ) J 



v n ,-k(r a ) 



(3) 



such that 



£u ni k(r a ) = -y^Jsu„,k(r a + 5) 

6 

- [fi + ftVa)] U nt k(r a ), 

Au„,fc(r Q ) = A (r a )u n:k (r a ) 

+ ^2 A s(r a )u n ^ k (r a + 6), 
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where the gap functions are defined by 

A (r Q ) = V Q (c^(r a )c L (r a )); 
As(r a ) = V\(c](r a + 5)c;(r a )). 



(4) 



(5) 



(G) 
(7) 



We have introduced the index a labeling the two basis 
points in the unit cell, the label k which is the Fourier 
inverse of R, and S = 0, f, f — R, —R connecting nearest 
neighbors with different basis indices. The equations (j^) 
are subject to the self-consistency requirements 



A (r Q ) = Vq M» |fc (r Q .)«* i _ fc (r a )taiih 



nk 



(8) 



A a (r a ) 



2 



[ u n,k{r a + S)v*_ k (r a ) 



iki) 



+ U„,k(r a )v* _ fe (r Q + 5) 



tanh 



f3e n „ 



(9) 



where the sum is over positive energy eigenvalues e n _k 
only. 

The orthorhombicity of YBCO is modeled by an 
anisotropy in the hopping parameters, reflecting the ip=. 
creased electronic mobility associated with the chains ;c3 
throughout the present work we use t v /t x — 1.5, which 
approximates the observed a-b anisotropy in the mag- 
netic penetration depth,E3 and Vq = —V\ = — 3t x . Twin 
boundaries of up to 4|f| width (corresponding to approx- 
imately 22A) are investigated, with t y = t x and /i 1 < 
within the twin. The largest system studied is 100|r| in 
length with 400 fc-states, or 80 000 sites; periodic bound- 
ary conditions are assumed throughout. 
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FIG. 1. The unit cell of the finite-size system for the BdG 
calculations is shown as a solid line superimposed on a square 
lattice. Long and short dashed lines represent twin bound- 
aries of width and |f|, respectively. Basis points are labeled 
by circles and squares. 

At a twin boundary of zero width and yJ = 0, we 
find within the BdG theory that for all temperatures the 
dominant rf-wave component of the order parameter is 
virtually unaffected. The extended and on-site s-wave 
components, whose bulk values are approximately 10% 
of A,j, go from their near-bulk values to zero, over a sin- 
gle lattice spacing r, reversing their sign relative to 
on either side of the boundary. As the impurity strength 
is increased at low temperatures, however, the d-wave 
and s-wave components become increasingly perturbed 
from their bulk values over the coherence length £d+ s (T), 
where £d+ s (0) ~ \f\ in the present work. When the mag- 
nitude of the d-wave component in the twin boundary 
is suppressed to approximately half its bulk value, an 
additional complex s-wave component may be nucleated 
near the twin edge, breaking time-reversal symmetry. We 
have found no evidence for a phase transition to a bulk 
T-violating state in a uniform system. 

The real and imaginary parts of the various compo- 
nents of the order parameter are shown in Fig. ^ for a 
twin boundary with fi 1 = — 10t x , fj, = —t x , T = 0, and 
boundary width Wt — 4|f|. While all the components 
go to zero rapidly within the twin boundary, both the 
real and imaginary parts of the s-wave gap functions are 
enhanced near the twin edge. In the immediate vicinity 
of the twin boundary, the real s-wave components are 
perturbed from their bulk values over a short distance 
comparable to ^+s(0), reflecting the local nucleation of 
additional s-wave components through spatial variations 
of the dominant d-wave component. The presence of fi- 
nite complex gap functions in the bulk, however, implies 
that the imaginary components vary over a different char- 
acteristic distance £j s (0) 3> £d +s (0). This longer length 
scale, as well as the comparable sizes of Im(Ao) near the 
twin edge and the bulk A^, is consistent with the GL 
predictionsQ discussed above. 
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FIG. 2. The (a) real and (b) complex components of the 
d-wave, on-site, and extended s-wave gap functions (corre- 
sponding to long-dashed, short-dashed, and solid lines respec- 
tively) are shown as a function of distance r perpendicular to 
the twin boundary. In (b), Im(A s ) is multiplied by a factor 
of —10 to facilitate comparison with Im(Ao). Results are ob- 
tained using fi = —t x , fj, 1 = — Wtx, T — 0, and twin width 
4|f|. Twin boundaries are centered at positions r — (which 
is equivalent to r = 100 by periodic boundary conditions) and 
r = 50 in units of If I. 
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FIG. 3. The maximum value of the complex on-site (dashed 
lines) and extended s-wave (solid lines) components relative 
to the bulk value of the d-wave order parameter are shown as a 
function of (a) \fj, | (note /j, < 0) and Wt at zero temperature 
and \i = — t x , and (b) T/T c and fi. In (a) thin and bold 
lines correspond to twin boundary widths VFr/|r | = and 
1, respectively. In (b), thin and bold lines denote results for 
fi = — 2t x and /i = — t x , and we have chosen fi 1 = — W0t x , 
and Wt = 0. 

As shown in Fig. [|, the size of the complex s-wave 
component nucleated near the twin boundary is ex- 
tremely sensitive to the temperature, impurity strength, 
and the width of the twin boundary. At zero temper- 



ature, Fig. 3(a), a T-violating state first appears for 
H 1 « —2.7t x at all electron densities for an impurity 
line (i.e. Wt — 0). As the impurity strength continues 
to increase, the perturbation of the <i-wave component, 
and the maximum values of the imaginary s-wave com- 
ponents, begin to saturate. For Wt > 0, however, a 
lower impurity strength can give rise to T-violation at 
zero temperature, since the d-wave component is already 
suppressed by approximately 20% in a locally tetragonal 
region (with t x = t y ) even for fi 1 = 0. Increasing Wt 
beyond approximately 3|f| has no further effect. This re- 
sult, valid for all electron densities, is also consistent with 
the GL prediction^ that local tetragonal symmetry could 
favor a time-reversal breaking state at low temperature. 

The growth of all the s-wave components with decreas- 
ing chemical potential reflects the impending instability 
of the system against bulk dominant s-wav£-.s«percon- 
ductivity at slightly lower electron densities Olla As the 
temperature is increased at finite fi 1 , the complex com- 
ponent decreases to zero as \J\ — T/T*; the transition 
temperature T* is strongly density-dependent, scaling 
roughly with A^. The same T* is obtained for wider 
twin boundaries at a given density, though the magni- 
tudes of the complex s-wave components increases with 
increasing Wt- 

The spatial variation of the s-wave component's phase 
relative to implied by Fig. || leads to currents which 
flow parallel to the twin surface and in opposite di- 
rections on either side of the twin boundary. The 
strong impurity potential therefore mimics a line of 
temperature-dependent magnetic flux passing through 
the twin boundary and oriented parallel to the c-axis. 
As shown by the differential conductance in Fig. 
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+ K/cMI 2 /' (v + E„ lfc )] 



(10) 



where /' is the voltage-derivative of a Fermi function, 
this effective field splits the low-energy band of virtual- 
bound states-associated with Andreev reflections at the 
twin surface.tll Alternatively, the presence of two low- 
energy quasiparticle peaks in the tunneling conductance 
at low temperatures reflects the existence of a physical 
gap in the excitation spectrum, proportional to the mag- 
nitude of the total complex s-wave component. The zero- 
temperature maximum peak-to-peak separation, found in 
Fig. [|a to be approximately 0.2t x ~ 2meV (where t x is 
estimated from T c = 0.51^2; rj 60 K at this carrier con- 
centration (n) w 0.75), grows with increasing \Vo\ and 
but diminishes with increasing temperature and dis- 
tance from the twin boundary. While the tunneling con- 
ductance exhibits low-temperature features that are no 
doubt finite-size effects, it is evident that the low energy 
band splits at a temperature 0.1T C « 6 K which is consid- 
erably lower than the T* ~ 0.5T C estimated from Fig. ^b. 
A comparable splitting of the zero-energy peak has been 



3 



recently observedQ in tunnelims-.spectra of YBCO sur- 
faces, and has been interpreted^ as a clear signature of 
T-violation. 
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FIG. 4. A low-energy portion of the tunneling conductance 
near a twin boundary is shown as a function of energy, dis- 
tance r from the twin edge where (a) through (c) correspond 
to |f | through 3|f while (d) illustrates the bulk, and tempera- 
tures T = (lower), T = 0.05T C (offset 0.1), T = 0.1T C (offset 
0.2) and T = 0.15T C (offset 0.3). Parameters are as in Fig. [| 

In summary, we have found evidence for time-reversal 
symmetry breaking near twin boundaries in a rf-wave 
orthorhombic superconductor at low temperatures T < 
T* < T c , where T*/T c scales approximately with the size 
of the bulk d-wave gap. The magnitudes of the complex 
s-wave components associated with the T-violation de- 
pend strongly on the chemical potential and depletion 
of the carrier density in the twin boundary. These s- 
wave gap functions could be responsible for the finite 
Josephson currents observed in c-axis tunnel junctions 
to heavily- twinned YBCO. til As a consequence of the 
time-reversal breaking, the low energy quasiparticle peak 
in the tunneling conductance (related to the zero-bias 
anomaly in STS) is predicted to split in the vicinity of 
the twin edge. 
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